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Abstract 

A bypass transition model has been implemented in the Wind-US Reynolds Averaged Navier- 
Stokes (RANS) solver. The model is based on the Shear Stress Transport (SST) turbulence 
model and was built starting from a previous SST-based transition model. Several modi- 
fications were made to enable (1) consistent solutions regardless of flow field initialization 
procedure and (2) fully turbulent flow beyond the transition region. This model is intended 
for flows where bypass transition, in which the transition process is dominated by large 
freestream disturbances, is the key transition mechanism as opposed to transition dictated 
by modal growth. Validation of the new transition model is performed for flows ranging 
from incompressible to hypersonic conditions. 


Nomenclature 

Cf skin friction coefficient 

k turbulent kinetic energy 

M Mach number 

M t turbulent Mach number 

P static pressure 

V k production of turbulent kinetic energy 

q w wall heat flux 

R t turbulent Reynolds number = pk/ pu 

R y turbulent Reynolds number based on wall distance = pyy/k/p 

Re x plate Reynolds number 

Reot Reynolds number based on momentum thickness 

Re u vorticity-based Reynolds number 

s streamline coordinate 

Sij rate of strain tensor 

St Stanton number 
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t time 

T 0 stagnation temperature 

Uj velocity 

x, y, z Cartesian coordinates 

y + wall normal coordinate 

k Von Karman constant 

H dynamic viscosity 

fit dynamic eddy viscosity 

v kinematic viscosity 

v t kinematic eddy viscosity 

u specific turbulent dissipation rate 

f] vorticity magnitude 

p density 

r turbulent time scale 

Tij viscous stress tensor 

tJa turbulent stress tensor 

oo freestream 


Introduction 

This report describes the implementation of an engineering transition model in Wind-US. 
Wind-US is the production Reynolds Averaged Navier-Stokes (RANS) solver of the NPARC 
Alliance, a formal partnership of NASA Glenn Research Center and the U.S. Air Force Arnold 
Engineering Development Center (AEDC), with significant participation by the Boeing Com- 
pany . 1 It is used by the NPARC Alliance partners and a number of other U.S. organizations 
spread across government, academia and industry to conduct analyses of a broad range of 
aerospace applications. As part of the ongoing development and improvement of the soft- 
ware for hypersonic simulations, as sponsored by the Office of the Secretary of Defense (OSD) 
Testing and Evaluation, Science and Technology (T&E/S&T) Program and the NASA Fun- 
damental Aeronautics Program, the desire for an accurate laminar-to-turbulent transition 
prediction capability was highlighted. 

The transition model development process began with a literature survey to select a base- 
line model that would fit the needs of engineers and scientists using the Wind-US solver. 
Once completed, a method of transition prediction was implemented and subsequently re- 
fined. After this initial code development and calibration stage, a number of experimental 
benchmarks were evaluated. This report describes the laminar-to-turbulent transition pre- 
diction methodology chosen and its justifications, the modifications necessary to calibrate 
the transition model using existing experimental data, and validation of the new model for 
benchmark transitional flows. 
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Motivation 


Laminar-to-turbulent transition is a field of continuing research. Despite decades of study 
there are still many aspects of the transition process that are not well understood. Given 
this, a universally applicable and sufficiently accurate method of determining the onset of 
transition is not possible for most flows of common engineering interest. This is in spite of 
the importance of the transition onset location to the calculation of drag, heat transfer, and 
other boundary layer characteristics. If the transition location is not well known, substantial 
error can be introduced into the simulation. Therefore, it is important in engineering and 
design applications to be able to predict transition location to a reasonable accuracy. 

The onset of transition is affected by numerous factors including the freestream flow 
properties such as turbulence level, Mach number, temperature, and pressure; as well as 
properties of the object itself such as surface temperature, surface finish, and its alignment 
relative to the flow. All of these factors make universal transition prediction a daunting 
challenge. While it is not currently possible to develop a scheme capable of taking all of 
these issues into account, it is possible to address specific areas and develop techniques for 
dealing with the instabilities caused by a more limited set of parameters. 

Transition occurs due to instabilities developed in the Navier-Stokes equations that govern 
the physics of fluid flow. A RANS solver like Wind-US is capable of mimicking the transition 
process (and some such as Wilcox, 2 argue that the growth of instabilities is replicated in a 
somewhat similar fashion) via a rapid growth in the production of turbulence at a point in the 
boundary layer. However, the most commonly used RANS techniques utilize a turbulence 
model solving one or two additional transport equations to represent the mean turbulent 
effect, and universally predict the onset of transition to occur much too close to the origin 
of the boundary layer. That is, they generate instabilities at an unrealistically low Reynolds 
number, and the transition to turbulence occurs well upstream of experimental data. The 
motivation of the present work is to implement a method that more accurately simulates the 
onset of this instability, moving the numerically predicted transition location downstream to 
a value consistent with experimental results. 

As explained above, a transition prediction scheme that is universally applicable is not 
currently possible. In order to narrow the scope of the problem, the present work will focus 
on what is known as “bypass” transition. Figure 1 shows a schematic “road map” of the 
transition process based on the work of Morkovin, Reshotko and Herbert. 3 On the left of 
the figure are the so-called natural transition mechanisms that arise from a linear stability 
analysis of the Navier-Stokes equations. As the magnitudes of the boundary layer distur- 
bances increase (moving to the right side of the figure), these linear modes are “bypassed” 
and transition occurs more rapidly due to non-linear interactions. Freestream turbulence 
is one of the mechanisms that generate these disturbances. The freestream turbulence in- 
tensity (FSTI) in the levels analyzed in the present work, and typical of large high-speed 
wind tunnel facilities, falls into this bypass category. The focus here will be on modifying 
the equations solved by Wind-US to accurately reproduce the behavior of flow undergoing 
bypass transition due to free stream turbulence. 

Bypass transition due to freestream turbulence has been studied with mixed success 
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for several decades. Large data sets evaluating the effect of increasing turbulence intensity 
on transition location have been compiled, though the agreement among them is mixed 4,5 
and their applicability is somewhat limited. In addition to the shortcomings of available 
incompressible data, very little reliable and systematic transition data correlating onset 
location with disturbance environment is available for supersonic and hypersonic flows. In 
many instances, only general trends derived from simple correlations can be found for these 
complicated cases. The goal of the current model development is to serve as an engineering 
approximation and the present work will attempt to devise a formulation that will reproduce 
a transition location that falls within the spread of available data. 


Turbulence Model 


The Wind-US CFD code contains several turbulence models for providing closure to the 
R.ANS equations. The work described in this paper utilizes one of the most commonly used 
turbulence models, the Shear Stress Transport (SST) model of Menter. 6 References 7 and 8 
have shown that the SST model generally provides the most accurate calculation of turbulent 
wall bounded flows of any of the one- or two-equation models in Wind-LIS. Equations 1-11 
show the formulation used by Ref. 9 for the SST model with a provision for simulating 
transition. The model provides closure to the RANS equations through the addition of two- 
turbulent transport equations for the turbulent kinetic energy, k, and the specific dissipation 
rate, u. This k — co model is active near viscous wall boundaries, and away from the wall 
the model transitions to a formulation based on a traditional k — e turbulence model. This 
is achieved by recasting the k — e based equations into a k — u formulation, and activating 
the additional terms by means of a switching function F\ . 
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and from these the turbulent viscosity is given as 


IM 

= min j 

( *P k aipk\ 

v" w : vf 2 J 

( 9 ) 

f 2 

= tanh(arg 2 ) 

(10) 

arg 2 

= min j 

f Vk 500zA 
i (3*ujy 1 uy 2 j 

(11) 


The constants (0 1 ) associated with these equations for the inner, k — u>, model are: 


<Jki = 0.85, 
CTwl = 0.5, 


0* = 0.09 


5/18 + (Rt[ 8) 4 

1 + (Rt/8) 


4 3 


ctl = 


.025 + Rt / 6 


0i = 0.075, 


1 1 + Rt / 6 ’ 

5 0.1 + R t /2.7 


OL i — 


9 1 + Rt/2.7 ’ 


a i = 0.31 


and the constants (0 2 ) for the transformed k — e model are: 

cr k 2 = 1.0, 02 = 0.09, al = 1, 

<7^2 = 0.856, 02 = 0.0828, a 2 = 0.4403, cq = 0.31 


These constants are blended using the same switching function, F\ that is found in the model 
equations such that 0 = i00i + (1 — 70)0 2 for any of the given parameters. 

The SST turbulence model has become popular for its ability to handle separated flows 
and complex geometry in the near wall region due to the strengths of the k — u) model, while 
maintaining the characteristics of the k — e model to be more accurate in free shear layers. 
Also, SST has shown to be more resistant to problems resulting from non-orthogonal grids 
that can result from complex geometries. 

The modification to these equations used to control the transition onset location is the 
“Production Term Modifier,” PTM, found in Eq. 1 which reduces the production of tur- 
bulent kinetic energy, Vk in the A:-equation. It is not used in the u equation. In order to 
properly model fully turbulent flow, PTA / I must be set to 1 in the fully turbulent regions. 
To prevent the transition to turbulence that occurs much too far upstream in the baseline 
SST model, PTM should be less than 1 in the laminar boundary layer. In this way, Vk 
is limited and the transition to full turbulence occurs less rapidly. An added stipulation is 
that the form of PTM only use local flow quantities and not be dependent on integrated 
parameters such as displacement or momentum thickness. Requiring integral quantities or 
other non-local values requires increased computational time and risks a loss of generality 
for complex flows with multiple viscous walls. By relying on local terms, these problems are 
avoided. 

The general form of PTAd is taken from the model initially presented in Ref. 9 and 
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modified slightly in Ref. 10. Equations 12-15 show the model as described in Ref. 10. 


PTM = 1 - 0.94(PTM1 + PTM2)F 3 


f 3 


= e 


( 


f) 


4 


PTM 1 


1 - [(3.28 x 10~ 4 )Re v - (3.94 x 10~ 7 ) Rel + (1.43 x KT 10 )^] ; 
1 - [0.12 + (1.00 x 10~ 5 )Pe„] ; 


PTM 2 



0-4 Re v . 
80 ’ 


K < 0 
K > 0 


( 12 ) 

(13) 

Re v < 1000 
Re v > 1000 

(14) 

(15) 


where the pressure gradient parameter, K, is given by 


K 


!± 

p 2 U 3 


'1 - M 2 ' 


dp 

ds 


(16) 


In these equations Rt is the turbulent Reynolds number and is effectively the ratio of 
turbulent to laminar viscosity, and Re v = py 2 Tl/ p is the vorticity based Reynolds number. 
Together, the claim in Refs. 9-10 is that these terms will limit Vk in the laminar region 
where it is over-predicted by SST and allow fully turbulent production downstream. 

Ref. 9 proposed a correlation in experimental data between the turbulent fluctuations 
upstream of transition and Re v . It was proposed to use this correlation as a means of 
controlling the location of Vk in the wall normal direction. The quantity Rt is used to 
control the PTM function in the flow direction as the freestream disturbances propagate 
into the boundary layer. The F 3 switching function acts as a gate, determining the threshold 
R t to begin turbulent production. The following sections address the performance of this 
initial model and the resultant changes that became necessary. In Refs. 9 and 10, PTM 2 
was formulated for flows with significant internal flow pressure gradients, and specifically for 
flows within low pressure turbine stages. In all cases examined in this work, where flows with 
significant pressure gradients were not examined, PTM 2 was not found to be significant, 
and as such no modifications to this parameter are addressed in the present work. 


Computational Setup 

For the initial transition model investigations, all computations were performed on a zero- 
pressure-gradient flat plate 200 in. (5.08 m) in length. This length was selected to give 
the appropriate range of transitional Reynolds numbers. A region containing 21 inviscid 
points was used upstream of the leading edge of the plate to ensure uniform flow conditions 
approaching the flat plate. The number of grid points was varied considerably to investigate 
grid convergence and will be discussed in later sections. The computational domain was 36 
in. (914.4 mm) high in the wall normal direction ending in a freestream condition. In all 
cases the axial points were clustered tightly in the transitional region, and the vertical points 
ensured a y + value of the first point off the wall of less than 1. 
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To baseline the performance of the proposed model, experimental data used by Refs. 
9 and 10 was evaluated. The data used for transition in the bypass regime is given by 
Ref. 11 for a flat plate with a turbulence intensity at the leading edge of 3.3% referred 
to as the T3A data. 11 The initial simulations evaluating the proposed model duplicated 
these non-dimensional experimental parameters as a starting point. However, the nature 
of the Wind-US solver makes it very difficult to exactly duplicate the T3A conditions as 
the data was taken at extremely low speeds which a strictly compressible flow solver, like 
Wind-US, does not handle easily. However the equivalent turbulence intensity and Reynolds 
numbers can be matched for a case where the speed is increased to a faster, but still nearly 
incompressible, Mach 0.2 freestream condition. 

Due to the increased freestream velocity, other parameters must be altered to give the 
appropriate non-dimensional conditions. To properly model the turbulence decay found in 
the T3A experiments, the input values of k and u need to be specified at the inlet. These two 
values describe the initial turbulence intensity as well its rate of decay, and are found, for the 
SST model, using equations 17 and 18, where the freestream (oo) quantities are the desired 
dimensional values at the leading edge needed to reproduce the non-dimensional FST1 found 
in Eq. 19, and the values of (3 and /3* are those of the outer model. 


v(s) = [Ps/Uoc + 1/uJoo] 1 (17) 

k(s) = /Coo [/^oos/Uoo + l]~ P //3 (18) 

FSTI( %) = ±Jf k (19) 


Note that s is the upstream axial position relative to the leading edge, and should have a 
negative value. These equations are used to ford the necessary upstream conditions for the 
simulation, ft may be observed that as one moves upstream, k and u get larger, and the 
dissipation rate equation (Eq. 17) in particular yields a maximum value for s that one can 
set an inflow upstream of the leading edge. Subsequent simulations at different FST1 levels 
were evaluated by changing k at the inlet. As discussed later in this report, the limiting 
behavior of the SST model shown in Eqs. 1 and 2 and the actual behavior of the SST model 
installed in Wind-US, replicate the decay of turbulence quite well. 


Convergence Behavior 

In the following section there will be frequent comparisons to the parameters used in the 
initial model. To facilitate easier comparison of the different parameters being investigated, 
the F 3 function defined in Eq. 13 will be referred to in the form shown in Eq. 20. 

F 3 = (20) 

Examining the baseline transition formulations, one Ends that a = 5 and (3 — 4 in the model 
in Ref. 10 and a = 6.5 in the original formulation of Ref. 9. The parameters a and f3 used 
here are not to be confused with the coefficients of the baseline SST model. 
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The transition onset location using the method as presented by Ref. 10 demonstrates 
considerable sensitivity to how the solution evolves and converges in computational time. 
For values of f3 > 2 and/or a > 3 the solution given by the transition model was found to 
depend strongly on whether the flow field was initialized from a converged turbulent flow, 
a converged laminar flow, or from uniform freestream conditions. Figure 2 shows one case 
where starting the solution from an initial laminar profile differed from the solution obtained 
by starting with an initially turbulent flow. Note that in this and subsequent figures, two 
correlation lines are provided on each skin friction plot, one for a fully laminar boundary 
layer and the other for a fully turbulent boundary layer, to compare with the transitional 
boundary layer solutions. These difficulties in obtaining consistent results regardless of the 
initialization procedure motivated efforts to modify the transition model. While Refs. 9 and 
10 attempt to address this issue by suggesting the transition model only be applied after a 
converged laminar solution is found, generalizing the model to be independent of start-up 
procedure was desirable for the implementation in Wind-US discussed here. 


Effect of Parameter a 

An analysis of the PT M 1 modifier in the PT M equation shows that it is relatively constant 
in the area of interest throughout the boundary layer in the pre-transition and transitional 
region, and plays a much smaller role than the F 3 switching function. Calibration of PTM 1 
will be discussed in subsequent sections. 

Varying the parameter a in Eq. 20 has a large effect on both the transition location and 
the behavior of the solution far downstream. The physical effect of a is to serve as a threshold 
for activating PTM. As Rt approaches a, PTM quickly switches from laminar (no turbulent 
production) to a state enabling turbulent production. The transition location is dictated by 
where this threshold is set. In addition to the undesirable convergence behavior discussed 
in the previous section, the suggested value of a = 5 from Ref. 10 not only locates the 
onset of transition significantly further downstream than predicted in experiments, but also 
causes the solution to never become fully turbulent (given this, a = 6.5 from Ref. 9 was not 
evaluated). The reason for this is that very near the wall in a fully turbulent boundary layer, 
R t is on the order of 5 or less into the beginning of the logarithmic region. Eq. 20 shows that 
the switching function for damping production, F 3l will prevent production of turbulence in 
the fully turbulent portion of the boundary layer, despite the fact that F 3 is only intended 
to be active prior to transition. Figure 3 shows the evolution of skin friction coefficient as 
the flow proceeds downstream. It may be observed that the skin friction coefficient for the 
flat-plate flow falls well short of the turbulent correlation 12 in the limit as Re x becomes large. 
Reducing the value of a to 3 improves the transition onset location somewhat. However, F 3 
is still damping the production of turbulence in the fully turbulent region, which results in 
a deficit in Cf compared to the expected value. 

The cause of this can be seen quantitatively by viewing contours of the Production 
Term Modifier (PTM) for the two simulations. Regardless of the behavior associated with 
predicting the onset of transition, PTM needs to become insignificant further downstream 
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so as not to influence the fully turbulent region. Figures 4(a) and 4(b) show that decreasing 
the value of a has a marked effect on how far PT M extends into the flow, but neither value 
fully eliminates the downstream influence. This is problematic in that the overwhelming 
majority of problems of interest will have fully turbulent flow over a significant portion of 
the calculation domain. Even if the model were perfectly calibrated with respect to the 
transition location, an error in the fully turbulent behavior would still result in significant 
error in the downstream region. It is essential that the model produce the correct behavior 
far downstream for it to be a usable tool. This is a fundamental failing of the baseline model 
that cannot be remedied with mere “tweaks” to the calibration, and a new approach must 
be developed. This new approach will be detailed in the following sections. 

Despite this shortcoming, and after significant numerical experimentation it was found 
that using a value of a = 3 produced the most consistent results in terms of the problems 
associated with the flow initialization. Starting from laminar, turbulent, or uniform flow 
fields all produced the same result using a = 3. This was an improvement over the baseline 
model formulation, and will be used as a starting point in further results. 


Effect of Parameter (3 

To analyze the effect of /3, a was held fixed at the value found in the previous section ( a = 3) 
and variations in / 3 were examined. The physical effect of f3 is to determine how rapidly PTM 
changes from the laminar to turbulent production mode as one deviates from the threshold 
value set by a. Increasing (3 makes that transition more abrupt, while decreasing (3 leads to 
a more gradual change. However, as the change becomes more gradual the flow experiences 
the effect of PTM over a larger range of R tl and lowering [3 too much produces inaccurate 
downstream behavior. (3 — 2 was found to be the most repeatable and numerically stable 
setting, as well as giving values in the neighborhood of the expected transition location. This 
result is shown in Figure 5. Although in this instance (3 — 2 and (3 — 4 are in agreement, 
this is not universally the case. (3 — 2 was used as a starting point for future work to make 
the transition process more gradual. 


Effect of Grid Resolution 

In order to study the effect of grid resolution, the optimized combination of a and (3 was 
chosen (a = 3, (3 — 2). Grids were created with two different numbers of points in the wall 
normal direction (j = 161 and 321) and a progressively increasing number of points in the 
stream- wise direction ( % = 230 through % = 1023). The present model was found to be quite 
sensitive to axial grid resolution, for reasons that are as yet unknown. Work by Rumsey 13 
documents similar sensitivity, and an evaluation of the unmodified SST model shows that its 
(much too early) transition location, in fully turbulent mode, is also strongly dependent on 
grid spacing although it is typically ignored. It should be noted that this is not a particular 
deficiency of the SST model as all two-equation models seem to demonstrate this behavior in 
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fully turbulent mode. Figure 6 shows the resolution needed to achieve grid independence. In 
the present work, over 700 axial points were found to be necessary for the baseline conditions 
of 3.5% FSTI and Mach 0.2. Only two values of wall normal spacing were evaluated as they 
produced consistent results. More specifically, the region well beyond the transition location 
(. Re x > 10 6 ) contained only 50 points, meaning over 650 points were required in the first 
initial part of the flat plate in order to achieve a grid independent answer. The points in 
the transition region were clustered such that Arc" 1 " = 62.8 calculated using a nominal Cf = 
0.003. 

For an engineering tool the grids required to obtain solutions independent of resolution 
and the associated computational time necessary are somewhat prohibitive. This sensitivity 
of transition location was observed by others (i.e. Ref. 13). In Ref. 14, the streamwise 
grid sensitivity was attributed to the first oder upwinding of the advection terms in the 
turbulence model. Our attempts to use higher order upwinding and TVD schemes were 
unsuccessful in reducing the grid sensitivity. The single biggest hurdle that remains in the 
use of this model will continue to be the density of the grid needed to accurately resolve the 
transition location. The results in the following sections should be viewed in context of this 
shortcoming. 


F3 Correction 

The preceding sections establish a baseline model that already exhibits substantial improve- 
ments over the initial form of Ref. 10. Specifically, it converges to a consistent solution 
regardless of initial values, and the transition onset is more consistent with experimental 
data. Once this base was established, further numerical experimentation was undertaken to 
improve overall accuracy and generality of the model. The sensitivity to grid spacing, current 
transition location, and the general shape of the transition region all pointed to the fact that 
a more gradual switching function, F :i . for controlling PTAd was necessary. However, this 
switching needs to be done without altering the correct limiting behavior already found for 
the F : > function. That is, F 3 should remain the same as Rt — > 0 and as Rt —■ ► 00 , but smooth 
the transition in between. In addition, the function should not remove the production of 
turbulence near the wall in the fully turbulent region. Two corrections to the form of F 3 to 
accomplish this are described in the subsequent sections. 

To moderate the behavior of F 3 in the transition region while maintaining its behavior 
at the extremes, the pass band of a Gauss filter was used to expand the values of R t that 
kept F 3 in an area of 50% production, while the stop band prevented any changes as R t — > 0 
and as R t — > 00 . The distribution was centered around a and expanded to produce the 
proper transition location. Table 1 shows the comparison between the initial model and the 
final value of F 3 that was found after numerical experimentation. This result is displayed 
graphically in Figure 7. The more gradual change in F 3 serves to decrease sensitivity and 
also to locate the transition onset properly. 

Questions may be raised regarding the deviation of the new F 3 modifier from the initial 
formulations. There is no reason to expect the shape to have any unforeseen effects on the 
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Table 1: F 3 Comparison 


behavior of solution despite this. Recall that the role of F 3 is to act as a gate, allowing the 
generation of turbulence above some threshold value of R t . The present functional form of 
F 3 simply increases the range over which moderate but not full production is taking place. 
Maintaining the exponential decay at the extremes prevents the function from introducing 
undesirable behavior. 


Fully Turbulent Correction 


The change in the prior section to the F 3 switching function is designed to remedy problems 
with the profile and transition location, but the downstream behavior in the fully turbulent 
region remains an issue. To prevent PTM from continuing to act in the downstream region, 
a limiter based on the non-dimensional coordinate y + is employed. 


y 



( 21 ) 


The parameter y + was chosen because it has predictable behavior relative to the turbulent 
boundary layer and is calculated on a point-wise basis, (i.e. it does not require any integral 
quantities to be calculated). The proposed downstream limiter is given by the form shown 
in Eq. 22 and was calibrated to eliminate the downstream remnants of PTM which are 
necessary in the laminar and transition regions but needs to be eliminated for the fully 
turbulent region. 

-^3 (modified) -^3 tanh ((2/+/17) 2 ) (22) 

Figure 8 shows the correction’s effectiveness in eliminating unwanted downstream behavior 
and allowing fully turbulent production. In these plots, the sharp decline in the contours 
indicates the transition location. Equation 22 introduces an additional dependence to the 
PTM formulation, which is now a function of R tl Re v and y + . 

One concern in using this form of limiter was that the formulation would not limit tur- 
bulent production in the laminar region near the wall where it might be necessary. However, 
evaluation of the flow held indicated that Vk is insignificant in this region when the how is 
laminar (due to low strain rate/vorticity, see Eq. 3) and thus was found to have no negative 
effects in the upstream laminar region. Downstream, it enables fully turbulent production 
and allows the skin friction prohles to accurately approach values specihed in various corre- 
lations. While relying on y + as a limiter deviates from the philosophy of the Menter SST 
turbulence model (and underlying Wilcox k — lo near wall model), the formulation enables 
the correct fully turbulent behavior and justifies its use here. 
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PTM 1 and Experimental Validation-Incompressible 

Given the extent to which the original model of Ref. 10 has been changed, the original 
formulation for PTM 1 was examined following the reformulation of F :i . The new modifier 
should still be based on Re v to exploit the properties outlined in Ref. 9 and properly account 
for freestream turbulence intensity. Specifically, a correlation between fluctuation intensity 
and Re v is selected to correctly locate Vk in the boundary layer. The modified expression 
that was investigated is as follows: 


PTM 1 


1 — Cptm 1 


[(3.28 x 10 -4 ).Re„ - (3.94 x 10~ 7 )Re 2 v + (1.43 x lO’ 10 )^] ; 
[0.12 + (1.00 x 10 -5 -Re„] ; 


Re v < 1000 
Re v > 1000 
(23) 


where Cptmi = 1.0 to yield the baseline formulation and the effect of increasing Cptmi was 
investigated as discussed next. 

To examine variation in Cptmi, we return to the T3A case described in previous sections. 
Unlike the simulations used for the benchmarks, where the T3A conditions served as an 
approximate starting point, the present simulation is matched to the precise turbulence 
intensity (3.3% at the leading edge) and decay rate as the experimental data. The SST 
model does a good job of reproducing the FSTI decay found in the T3A wind tunnel as 
the flow progresses downstream. Figure 9 shows a comparison between the decay found in 
Ref. 11 and the decay in the simulation. Despite the fairly steep gradients found here, the 
decay rate and initial value are successfully reproduced. Note that this plot shows a slight 
undershoot on the freestream intensity, which leads to an expectation that the transition 
onset in the simulation will occur slightly downstream of where it occurs in the experiment. 

Figure 10 shows a comparison of Cf for simulations of the T3A flat plate with and without 
the transition model. The fully turbulent SST solution shows the more rapid transition 
(much closer to a fully turbulent flow) that provided the impetus for the present work. It 
does a poor job of accurately capturing the early laminar behavior and the transition location 
which would lead to incorrect results for the drag on this flat plate. For the transition model 
solutions, the effect of Cptmi is to control the transition onset location. As with nearly 
all RANS-based transition prediction schemes, the width of the calculated transition region 
is shorter than indicated by experimental data. The solution with Cptmi = 2.0 captures 
the transition onset location best while decreasing Cptmi to 1.0 (default value from original 
formulation) captures the end of the transition zone (i.e. where fully turbulent flow is 
realized) best. 

Figure 11 shows the results of simulations for the T3A test case where several freestream 
intensities are investigated using Cptmi = 1-0 and Cptmi = 2.0, the bounds on the rec- 
ommended values of Cptmi- bi this figure, the ordinate is the transition Reynolds number 
based on 6, the momentum thickness. The plot is presented in this manner because the most 
frequently used engineering correlations are calibrated based on this method. In Fig. 11 the 
transition onset location is defined as the minimum C/, and the FSTI is defined midway 
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between the leading edge and the transition location to keep in convention with how the 
correlations were formulated. 4,5 This method of reporting the FST1 is an attempt to account 
for the decay of the turbulence in the tunnel in a simple and useful way. To maintain con- 
sistency, the numerical results are also plotted in this way. The results shown in Fig. 11 
indicate that the Cptmi — 2.0 solutions match the onset of transition location somewhat 
better that Cptmi = 1-0. Recall that the results of Fig. 10 showed that Cptmi = 1-0 
matches the end of transition region (fully turbulent) location better. 

The range of applicability demonstrated here makes the model appropriate within the 
regime of “bypass” transition flow situations introduced earlier. While the model will likely 
have difficulty handling freestream intensities that approach zero, as well as very high distur- 
bance environments, it should perform well in predictions for a range of conventional wind 
tunnels. One comment about the freestream intensity is necessary here. The correlations 
used as experimental validation for the model are based upon grid-generated turbulence. At 
a point in the wind tunnel upstream of the leading edge of the plate a disturbance gener- 
ator is placed in the flow to create isotropic turbulence whose magnitude can be properly 
controlled. Because these are disturbances placed into an otherwise “quiet” wind-tunnel, 
the FSTI decays in the stream-wise direction. In contrast, turbulence found in conventional 
wind-tunnels comes from a variety of sources including acoustics and inflow irregularities 
which may not decay as the flow proceeds down the tunnel. While the transition model 
itself would likely have little problem handling a constant FSTI as opposed to the decaying 
FSTI values used in the current simulations, the limiting form of the model has difficulty 
replicating this constant-turbulence behavior. More work is necessary to address the FSTI 
decay rate in the tunnels being simulated, and if this value is a constant, additional terms 
may be necessary in the turbulence model to generate this disturbance at an appropriate 
level. 

In concluding this section, we summarize the final form of the transition model that is 
recommended based on the modifications discussed in this report: 


PTM 

F 3 

P(Rt ) 

PTMl 


1 - 0.94(PTM1 + PTM 2) F 3 tanh ((y + / 17) 2 ) 

e-WRl - P(R,)) C-PiR,) 

2.5 ~ ( fit— 3 ) 2 

P 2 


[(3.28 E - 4 )Re v - (3.94A - 7 )Re 2 v + (1.43 E - 10)7^] ; 
[0.12 + (IE - 5 )Re v \ ; 


1.0 < Cptmi < 2.0 ( Recommended range ) 


(24) 

(25) 

(26) 

Re v < 1000 
Re v > 1000 

(27) 

(28) 


In Wind-US, the default value for Cptmi is set to 1. 
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Experimental Validation-Compressible 

Experimental data systematically analyzing the effect of freestream turbulence in higher 
Mach number flows is not as prevalent as in the incompressible case. Concrete correlations 
for higher Mach numbers equivalent to those used in the incompressible cases are not readily 
available. The most frequently used tools in current use involve single point correlations of 
the form Reg/M = constant. These correlations have been developed in free flight conditions 
and in conventional supersonic tunnels and can be used to establish trends regarding the 
transition behavior at varying Mach numbers. The accuracy of these correlations is uncer- 
tain and comparing the model to them in any quantitative way is difficult, but qualitative 
comparisons are useful and are examined here. 

A series of cases was examined for Mach 2.5, 3.5, and 4.5 to investigate the performance 
of the transition model in supersonic flow. A computational grid having 661 axial points 
and 161 points in the wall normal direction was used. These supersonic cases demonstrated 
lower grid sensitivity than for the incompressible cases in that solutions obtained with every 
other grid point in each direction were very similar to those using the full grid. Since we 
are interested in comparing locations of transition onset, the simulations used Cptmi = 2.0. 
Figure 12 shows a plot similar to that used for the incompressible calibration, with the lines 
of Reg/M = 100 plotted for the Mach numbers being simulated. 15 The shaded region shows 
the same expected transition band from the incompressible correlations of Fig. 11 only 
for reference. The results for these supersonic cases show that in the range of turbulence 
intensities common in supersonic wind-tunnels, the model is able to qualitatively reproduce 
the experimentally determined behavior. The parameters for the numerical simulation are 
the same as in the previous section, changing only the Mach number and initial turbulent 
intensity. The model is able to successfully reproduce the known stabilizing influence of 
higher Mach number flows. 

While none of the numerical data points correspond precisely to the correlations shown 
in Fig. 12, it has been documented that the model varies smoothly with increasing and 
decreasing intensities. It would not be difficult, nor would it be any more revealing, to find the 
exact FSTI that reproduces the transition location found with these simple correlations. In 
addition, this method of correlating the transition location has been under scrutiny. 16 These 
correlations are used here merely to establish a general trend and the relative magnitude of 
the expected changes due to compressibility effects. The lack of other meaningful empirical 
comparisons makes more sophisticated analysis difficult. 


Experimental Validation-Hypersonic 

Hypersonic flows present unique challenges in both modeling and experimentation. To base- 
line the present model’s ability to predict transitional behavior in hypersonic conditions, it 
was validated against transition data taken in the AEDC tunnel B on sharp nose cones at 
Mach 7.93. 17 The simulation is performed on a 7° half angle cone, 40 inches (1016.3 mm) in 
length, in agreement with Ref. 17, and as in the experiment, the wall temperature was set at 
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0.42X1), where X 0 is the freestream stagnation temperature. Several unit Reynolds numbers 
(Re/m) were evaluated in Ref. 17 to provide a complete scan of the transition region. These 
data show good agreement across a range of Re/m , allowing the simulation to use only one, 
6.8 x 10 6 Re/m, corresponding closely with the center of the experimental range. As in 
previous cases, the grid was clustered near the transition location and used 338 axial points 
and 161 points in the wall normal direction, maintaining a y + of the first point off the wall 
less than 1. Inflow conditions are calculated to match those of Ref. 17 and insight by Ref. 
18 that an inlet FST1 of 1.25% was appropriate for AEDC tunnel B. 

Figure 13 shows static pressure contours to validate the development of the shock around 
the cone. The shock angle is in very close agreement with theoretical predictions and demon- 
strates there is no need for further grid refinement to accurately capture the shock behavior. 
Figure 14 shows the Stanton number (St) as defined by Ref. 17 (and as shown in Eq. 29) ver- 
sus Reynolds number for the SST model alone, the SST transition model using Cptmi = 1-0 
and Cptmi = 2.0, and the experimental data. 

St = q w / (pooUoo(h(T 0 ) - h(T w )) (29) 

Although as in the incompressible case the behavior in the transition region itself is overly 
abrupt, the model accurately captures the minimum heat transfer value, corresponding to 
transition onset for Cptmi = 2.0 and the location where the follow becomes fully turbulent 
for Cptmi = 1-0. 

Of additional interest is the agreement of the simulation with the pre-transition data. 
The SST transition model accurately captures the laminar heat transfer behavior in the 
region leading up to transition. This is important in demonstrating the utility of the model, 
as finding a correct transition location would be irrelevant if the model was unable to produce 
the accurate laminar behavior prior to onset. Additionally, this is a non-trivial validation as 
the y + limiter placed in the PT M function allows for some turbulent production in the near 
wall region, even when the flow is laminar. However, the agreement between the transition 
model and the data verifies that this modifier has no noticeable effect on the flow behavior 
in the laminar region, even in the hypersonic case. Although no specific experimental case 
is shown for supersonic validation, the ability of the model to match behavior at subsonic 
and hypersonic Mach numbers demonstrates consistency across different flow regimes. 

This is a marked improvement over the SST model alone, which vastly over-predicts the 
heat transfer by indicating a fully turbulent state from the leading edge of the test article. 
For a model in a high-speed wind-tunnel that has a significant laminar region, the error due 
to transition onset can be quite large. The total heat transfer integrated over the length 
of the cone differs by 38.7% between the SST model and the transitional SST model with 
Cptmi = 1.0. Accurately capturing this behavior is especially important when trying to 
evaluate thermal properties and heat transfer behavior near the tip of a hypersonic vehicle. 
The large discrepancy in heat transfer rates indicated by the two models would substantially 
alter the vehicle’s predicted temperature profile. 
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Conclusions 


A laminar-to-turbulent transition method has been implemented and validated in the Wind- 
US flow solver. This method was built starting from a previous formulation centered around 
the Menter SST turbulence model, with significant modifications enabling a more accurate 
engineering calculation tool. The transition prediction model in the present work demon- 
strates several improvements over the previous form presented in Ref. 10. The current form 
of the “Production Term Modifier” prevents problematic convergence behavior, allows fully 
turbulent production downstream of the transition location, and is found to reproduce tran- 
sition onset accurately in a wide range of zero or slight pressure gradient flows. The model 
is intended for flows where bypass transition is the dominant laminar-to-turbulent transition 
mechanism, and is validated against experimental data for cases ranging from incompressible 
to hypersonic over a range of inlet turbulence intensities. While the model remains overly 
grid sensitive, and does not exactly replicate the behavior in the transitional region between 
fully laminar and fully turbulent, it nevertheless constitutes an improved engineering tool 
for finding behavior in flows where transitional effects are significant. Further, it accom- 
plishes this without the expense of adding any additional transport equations or requiring 
calculation of integrated parameters such as momentum thickness. This allows the model 
to be more efficient and avoid the ambiguities inherent when integrating in the wall normal 
direction with complex geometries. 

This report serves as a starting point for incorporating more accurate transition prediction 
into the Wind-US solver. More work into models involving additional transport equations, 
pressure gradients, and natural transition effects is necessary to cover the broad scope of 
transition phenomena. Additionally, efforts to eliminate the grid sensitivity that remains 
in the present model should be explored. Despite these limitations, this transition model 
represents an improvement in the ability of computational methods to complement wind- 
tunnel testing at a range of Mach numbers. 
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1 Figures 


Increasing Disturbance Level 



Figure 1: Transition road map, from Ref. 3 
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Figure 2: Spurious solutions demonstrating initialization effect 
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Figure 3: Dependence of transition onset on a 



Figure 4: PTM Contours 
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Figure 5: Dependence of transition onset on f3 
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Figure 6: Grid convergence 
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Figure 7: Comparison of original and modified F :i functional forms 



Figure 8: Comparison of F :i contours with and without y+ limiter 
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Figure 9: FSTI comparison 
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Figure 12: Supersonic Mach number comparison 
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